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Abstract: Type OB string theory has been proposed as the dual description of 
non-supersymmetric SU(N) Yang-Mills theory coupled to six scalars, in four dimen- 
sions. We study numerically and analytically the equations of motion of type OB 
gravity and we find RG trajectories of the dual theory that flow from an asymptot- 
ically free UV regime to a confining IR regime. In the UV we find a one-parameter 
family of solutions that approach asymptotically AdS^ x S 5 with a logarithmic flow 
of the coupling plus non-perturbative terms that correctly reproduce all UV and IR 
renormalon singularities. The first UV renormalon gives a contribution ~ Fi(E)/E 2 
and we are able to predict also the form of the function F\{E\ which, from the YM 
^ ■ side, corresponds to summing all multiple-chain bubble graphs. The fact that the 

positions of the renormalon singularities in the Borel plane come out correctly is a 
non-trivial test of the conjectured duality. 
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1. Introduction and summary 

The attempts to generalize the Maldacena conjecture |], 0, ||, |j to four- dimensional 
non-supersymmetric YM theories have followed different directions. The approach 
proposed by Witten is to start from the duality between type II strings and a five- 
dimensional super- YM; compactifying the fifth dimension on a circle with different 
boundary conditions for bosons and fermions one breaks supersymmetry and is left 
with an effective four- dimensional non-supersymmetric theory [[5]. A conceptually 
similar possibility is to reduce the supersymmetry to N = 1 or even N = with the 
addition of mass deformations || [7|. §. 

A different strategy is to start directly from a string theory with worldsheet 
supersymmetry but no spacetime supersymmetry. Polyakov has advocated this ap- 
proach and has suggested the use of non-critical strings 0, |I0 |, and Klebanov and 



Tseytlin [[11], [12| have adapted this suggestion to critical type OB string theory. On 
the gravity side one considers a configuration of N electric D3-brane in the type OB 
theory, and the dual field theory is conjectured to be SU(N) four- dimensional YM 
plus six adjoint scalars. The six scalars are the price to pay for working in the critical 
dimension. 

In this paper we follow the latter approach and we study the renormalization 
group (RG) flow derived from the equations of motion of type OB gravity. Klebanov 
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and Tseytlin |12[ have found that in the UV one recovers asymptotic freedom, 



IE E 

— — — ~ In const. In In — . (1.1) 

9ym{E) A A 

The numerical value of the beta function coefficients are not under control because 
of various corrections discussed below, but it is very encouraging to see the correct 



logarithmic behaviour. Minahan |13|, [14| has instead presented confining IR solutions. 



Further related work appeared in refs. [15|-|2~0fl. It is therefore natural to ask whether 



asymptotic freedom and confinement are connected by a single RG flow. 

In this paper we perform a detailed analytic and numerical investigation of the 
RG flow, and we find the following results. 



The UV solution found in refs. [12[ |1J] is just a member of a one-parameter 



family of solutions. The free parameter k, which is related to subleading terms 
in the UV, determines the IR behaviour of the solution. For k smaller than a 
critical value k c (that we can determine analytically) the solution flows in the 
IR toward a confining solution. The case considered in |l2j corresponds instead 



to k > k c , and it was found in [12| that in this case the solution flows in the 



IR toward a conformal fixed point at infinite coupling. 

We will then show, both numerically and analytically, the existence of RG 
trajectories that connect asymptotic freedom in the UV with confinement in 
the IR. 

We will find that in the UV, beside the logarithmic terms shown in eq. ( |1 . 1| ) , 
type OB gravity predicts also power corrections, 

1 , E , , E ^ /r ^A 2 ^,^A 4 A 4 1 . , . 
\n--c 1 \n\n- + c 2 F 1 (E)— + c 3 F 2 (E)— + c 4 — + 0(—) (1.2) 



g\ M {E) A A ' E 2 * zv 7 £ 4 4 £ 4 v £ 6 ' 

where the q are constants to be defined below. These power corrections match 
exactly the non-perturbative contributions coming from renormalons. The term 
~ Fi(E)/E 2 on the field theory side comes from the first UV renormalon, 
the term ~ F 2 {E)/E A from the second UV renormalon, and similarly we get 
the contributions of all UV renormalons. The term c^/E 4 is the leading IR 
renormalon, and we correctly find no IR renormalon contribution 1/E 2 . Thus, 
we get the right positions of all renormalon singularities in the Borel plane. 

Even the energy dependence of the prefactors, that is Fi(E) for the first UV 
renormalon and a constant for the leading IR renormalon, matches what is 
known from the field theory side. In particular, we will predict the energy 
dependence 

F 1 (E) ~ sin (Pt^-j^ + a T ^j (1.3) 
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with cut, Pt constants. Expanding at first order in (3 T we get a result consis- 
tent with what is obtained in QCD summing diagrams with a single bubble 
chain. The above closed form is a prediction of the dual string theory for the 
resummation of multiple chain bubble graphs, which would be very interesting 
to check from the field theory side. 

The above results hold at lowest order in a' and string loop corrections, and are 
presented in sects. 3-6. The fact that some features of the UV limit of the YM theory 
can be obtained from the gravity side is surprising, since the SU(N) YM and the 
gravity descriptions are duals to each other and are based on opposite expansions, 
in g 2 N and 1 / \fg 2 N , respectively. It is therefore important to study the corrections 
to the gravity predictions in the UV limit. We will discuss them in sect. 7. It has 
been found in ref. |12| that, if the a' corrections can be expressed in terms of the 
Weyl tensor, then the only price to pay in the UV is that the numerical values of the 
beta function coefficients cannot be predicted, but the logarithmic running of the 
coupling is unaltered. Furthermore we will find that, under the same conditions, the 
position of the renormalon singularities in the Borel plane is independent of the a' 
corrections. 



2. General strategy for computing the RG flow 

In this section we discuss the general strategy and the assumptions used to extract 
the RG flow in the dual theory from the gravity equations of motion of type strings. 
The low energy effective action of type OB strings pi , 22] , |23| is written in terms of 
a (closed string) tachyon field T, the dilaton 0, the metric, and a doubled set of RR 
fields. The terms involving the tachyon in the effective action have been discussed 
by Klebanov and Tseytlin [pi]] , and are 

J d 10 xVGe- 2 <f> (^G mn d m Td n T + V(T)^j + J d 10 xVG^f(T)\F 5 \ 2 . (2.1) 

Because of world-sheet supersymmetry, the tachyon potential V(T) has only even 
powers of T. We write (to have the same notations of [[H]) V(T) = —2g(T), and 

g (T) = - L -m 2 T 2 - AT 4 + 0(T 6 ) = ^T 2 + 0(T 4 ) . (2.2) 

The tachyon mass is m 2 = —2/a' and we have set a' = 1. The function f(T) 
describes the coupling of the tachyon to the 5- form field strength F 5 , and is given by 

f{T) = 1 + T+ -T 2 + 0(T 3 ) (2.3) 
2 



and the coupling to the other RR fields can be found in ref. (TT|. Note that, due to 
the doubling of RR fields, the 5-form field strength is not constrained to be self-dual, 
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so that I-F5I is non vanishing. D-branes in this theory have been studied in |24|, [II 
Another consequence of the doubling of the RR fields is that there are two kinds of 
D-branes, electrically and magnetically charged. It is then interesting to consider 
the solution of the gravity equations of motion corresponding to a large number N of 
coinciding electrically charged 3-branes, since from eq. ( [2.1| ) we see that the coupling 
of the 3-brane to the 5-form RR field strength can stabilize the closed string tachyon. 
The tachyon from open strings which ends on parallel D-branes is instead eliminated 



ds% = e^dp 2 + e~Ux 2 (4) + e^dttl*, . (2.4) 



by the GSO projection |TT], ||. In the Einstein frame, the ansatz for the 3-brane 
metric can be conveniently written in terms of two functions £ and rj of a transverse 
variable p, 

For the 4- form RR field Cu\ one takes Coi23(p) as the only non- vanishing com- 
ponent, so that F 5 = cfC(4) has only the component (-F 5 ) i23p; similarly one takes 
T = T(p),0 = 0(p). The equation of motion for the RR field can be explicitly 
integrated and gives an integration constant Q which is just the RR charge, and 
is proportional to N. The remaining equations of motion then separate into four 



dynamical equations [11, 14 



i + l -g{T)e^~^ = , (2.5) 

e + ^(T)e^ + |^e^ = 0, (2.6) 

77 + \g(T)e*P-*" + 8e- 4?? = , (2.7) 

f + ^PP^ + 2</(T)e^- 5 " = , (2.8) 

and (due to the invariance under reparametrizations of p) a constraint on the initial 
values, conserved by the dynamical equations, 

f + I - 5r) 2 + ^ + 9(T)e^ + 20e- 4 " - Q 2 j^ = . (2.9) 

The dot is the derivative with respect to p. The explicit form of f(T) plays an 
important role in the UV regime, where the tachyon is stabilized at its minimum, 
/'(To) = Q. 

There are uncertainties due to the fact that neither f(T) nor g(T) are known 
in closed form; actually the latter is not even unambiguously defined, since, using 
the equations of motions, a term ~ V 2 T in the effective action can be traded for a 
term ~ m 2 T so that, for instance, a higher derivative term ~ (V 2 T) 4 in the effective 
action can be reabsorbed into a ~ T 4 term. In the following we will at first set 
g(T) = (1/2)T 2 and f(T) = 1 + T + T 2 /2 (the minimum of f(T) is in this case at 
T = —1). We will then discuss how our results depend on these choices. 
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Our aim is to solve eqs. (|2.5|) - (|2.9|) , and to translate the solutions into a renor- 
malization group flow of the coupling g\ M °f the dual theory. The first issue is how 



to read g Y M from the gravity side. We follow the approach of |L3], [LJ, [L2[| and we 



read it from the quark-antiquark potential, computed from the Wilson loops as in 



refs. [26]]. This means that we take 

g 2 YM ~ exp{0(p)/2} (2.10) 
as our definition of the coupling. The exact proportionality factor will not concern us 



here since, as found in refs. [[14], [12] and as we will discuss below, the gravity solution 
is in any case subject to a number of uncertainties, due to a' corrections and to 
the exact form of f(T),g(T), that affect the numerical value of the proportionality 
coefficient in eq. fl2.10Q . 



The definition ( |2.10|) is apparently in contradiction with what one would obtain 



from the D-brane effective action, which seems to give g Y u ~ ex P{0} rather than 



exp{0/2}. However, as discussed in [fL2| , in the type theory there exists also 



a tachyon tadpole on the D-brane, so that the effective action of the D-brane is 
proportional not simply to but rather to k(T)e~^, with k(T) = 1 + T/4 + 0(T 2 ). 
Since in the UV the tachyon runs toward its minimum at T = —1, higher powers in 
k{T) cannot be neglected and the closed form of k{T) is needed to draw conclusions. 
It is clear, however, that this is a point which deserves further investigations. 

Eq. (|2.10| ) provides us with the dependence of the coupling on the transverse 
coordinate p. To obtain an RG flow, the second step is to connect p with a physical 



energy scale. This can be done as follows. From eq. ( |2.4j ) we see that a dilatation of 
the 4-dimensional coordinates, x — > Ax, can be reabsorbed into a rescaling of £ such 
that exp{— £/2} scales like an energy squared (while rescaling r\ and p so to keep the 
other terms invariant). We will then define the energy scale from exp{— £(p)/2} ~ 
E 2 , or 

In | = -\t{p)- (2.H) 

A is a scale that we will determine later. In the UV region the solutions will ap- 
proach AdSs x S 5 and, writing the metric of AdS$ as ~ du 2 /u 2 + v?dx 2 ,^ we see 



that exp(— £/2) ~ u and then the prescription ( |2.11| ) reduces to the by now stan- 
dard identification u ~ E. Far from the UV there is a certain arbitrariness in this 
definition, that we regard as corresponding to the arbitrariness in the choice of the 
renormalization scheme on the gauge theory side. Of course the details of the beta 
function are not independent of these choices, and only universal properties, like the 
existence of zeros and how the zeros are approached are what really matters. 

Solving the equations of motion we get (p(p) and £(p). Combining eqs. ( |2.10| ) 
and ( [2. 11| ) then provides the dependence of g\ M from the energy scale and therefore 
the RG flow. 
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3. The UV and IR asymptotics 



An asymptotic solution valid in the UV limit has been presented in ref. and, as 
a more systematic expansion, in ref. [jEJ . Direct inspection of eqs. ( f2.5| )-( pi?P reveals 
however that there is actually a one-parameter family of solutions. It is convenient 
to introduce a new variable y from 

p = e- y . (3.1) 

In the limit y — > oo the solution is 

1 /ln 2 u\ 
= -2 toy + 151n2 + -(39\ny + k) + — / , (3.2) 

V \ V J 

£ = -y + ^ 2 + - + -L(39 hxy + k - 104) + O (^f) , (3.3) 

y 2y V v J 

V 11 f\n 2 y\ 

r) = -f + ln2 + - + — 391ny + k - 38) + O l—f- , (3.4) 

2 y 2y 2 \ V J 

8 4 /ln 2 y\ 
T = -l + - + — (391n 2 / + /t-20) + O — y - . (3.5) 

y r V v J 

Here and in the following we set Q — 1; we see from the equations that the solution 
for generic Q can be recovered from 

=f(y;Q = l)+lnQ, 
<P(y;Q) = <P(y;Q = l)-hiQ, (3.6) 

while r^, T are independent of Q; k is a free parameter, and the solution found in 
ref. Jl2| corresponds to k = 0. In spite of the fact that in the UV k only appears in 



terms which look quite subleading, we will find that its value is very important for 
determining the IR behaviour of the solution. 1 In the large y limit £ ~ —y — > — oo. 
Eq. ( |2.11| ) then shows that y —>■ oo (or p — > 0) is the UV region. At y = oo the 
metric reduces to AdS 5 x S 5 . 

Using the prescription discussed in sect. 2, one can now extract the RG flow in 



the UV, with the result 12 



^ m(E) ~ ln(E/A)-f lnln(E/A) ' (3 ' 7) 



^^This free parameter comes out because, inserting into the equations of motion an ansatz of 



the form of eqs. (3.2)-( |3.5| ) with generic coefficients, one finds that in eq. (2.5) the terms of order 
(hiy)/y 2 cancel automatically, and therefore impose no constraint on the coefficients of the solution. 
Looking only to the terms up to 0(1/ y) one might think that this parameter can be removed by 
a conformal rescaling, but this is not true anymore when one considers also the terms 0(1/ y 2 ) in 
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Note that the logarithmic terms are independent of k. Eq. ( |3.7| ) should be compared 
to the running of the coupling in the proposed dual theory, that is four- dimensional 
SU (N) YM theory with 6 scalars in the adjoint representation, which at the two-loop 
level is 

9ym( e ) = —, —7 r- (3.8) 

b 1 (\n(E/A) + ^\n\n(E/A)) 



with 

We see that eq. ( |3.7| ) qualitatively reproduces the logarithmic running of the coupling 
in the UV. The precise value of the beta function coefficients cannot be reliably 
estimated, since the gravity prediction is affected by a' corrections (see refs. fl2], |H1 
and sect. [7|). Furthermore, with a generic tachyon potential g(T), and a generic 
function f(T) with a minimum at T = T , the factor 39/8 in eq. ( |3.7| ) is modified 
as |L4 

s i^i+(?m) 2 (3 io) 

so it is clear that the values of the numerical coefficients are not under control. 

In the opposite limit, p — > 00, it is again possible to find asymptotic solutions. 
In this case however there is a much larger variety of solutions. In particular Mina- 



han 14] has considered the generic behaviour 



<j> ~ (f)ip + 0o 

£^£iP + £o (3.11) 
77 ~ rjip + r]o 
T ~ tip + t , 



At large p this is a solution of the equations of motion if 

&,7&>0, 5»7i - - > (3.12) 



\<t>\ + \ei-H + \t\ = U- (3-13) 



and 



The latter equation follows from the constraint equation fl2.9p , while the inequalities 
( |3.12| ) ensure that all exponentials in eqs. (|2.5|)-(p^) are suppressed and therefore 
that the equations of motions are satisfied. 

In the limit p —>■ 00 we now have £ — > +00 and eq. ( |2 . 1 1| ) tells that this is the 
IR limit. To investigate confinement one can compute the quark-antiquark potential 
as in ref. p5|, E6l E7L EHJ. One considers a Wilson loop on the boundary p = 0, with 
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edges along the directions x, t, and looks for the classical string world-sheet p(x, t) 
which has the Wilson loop as its boundary. The Nambu-Goto action is 

1 



S 



NG 



2na' 



dadrJdet{G M Nd a X M d p X N ) 



(3.14) 



where Gmn is the string frame metric, which in ten dimensions is related to the 
Einstein frame metric Gmn by Gmn = ^ 2 Gmni anc ^ ^m! f° r the 3-brane solution 
is given in eq. ( |2.4|) . Inserting Gmn into eq. ( |3.14|) one finds that the static potential 
V(L) of a qq pair separated by a distance L along the direction x is obtained from |L4] 

1/2 



V(L) 



1 



2na' 



dx 



3 0— 5?7 




(3.15) 



after subtracting the divergent part corresponding to the energy of two separated 
massive quarks [^5|, £|. The above result holds for generic Q. The Q dependence can 
be extracted using eq. fl3.6|) and gives 



V(L) 



2na'y/Q Jo 



dx 




1/2 



The result depends crucially on the function 



; 2 (p) 



\e v ? Q= i — ► e l 



-£i)p+W>o-£o) 



(3.16) 



(3.17) 



(We define O , £o as the coefficient of the solution with Q — 1). If <pi < £x, / vanishes 
at p = oo, and the classical string configuration will go all the way to the region 
where p — >• oo, where there is no cost in energy in separating the quarks further; 
for 0i = £i, the classical string configuration again goes to the minimum of f(p) at 
p — > oo, but now /(oo) = exp{(0 o _ ^o)} is non- vanishing; the qq potential at large 
L then becomes 

V(L)~aL } (3.18) 

with 

a = - _ e ^-«°)/ 2 (3.19) 



If instead 0i > ^i, then the function f(p) diverges at p — > oo. On the other hand, it 
also diverges in the UV, where (0 — — > y — > +oo, see eqs. ( p.2|) -( p73|) . This means 
that it has a minimum at some value p m i n , where the classical string configuration 
goes. Then we have again a linear qq potential with 

1 



a 



f(pmin,Q = 1) 



(3.20) 



2-KOi'Q 

This gives the relation between the 'QCD' string tension a, the type string tension 
l/(27ra/) and the RR charge Q. In the case of SU(3) YM, the experiment gives 
a ~ (1 fm)~ 2 and eqs. ( |3.19| ) or ( |3.20| ) fix the value of a' for the dual string theory. 
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Summarizing, the condition for confinement is 



0i > 6 • (3.21) 

A very different 1R solution has been found by Klebanov and Tseytlin . In this 
case, approaching the IR limit, the Einstein frame metric becomes again asymptotic 
to AdS$ x S 5 , while the dilaton diverges and the tachyon goes to zero. Therefore the 
dual theory flows toward a conformally invariant point with infinite coupling. 

It is at this point natural to ask whether the one-parameter family of UV solu- 
tions discussed before is smoothly connected to any of these IR solutions. This is the 
issue that we will address in the next section numerically and in sect. [5] analytically. 



4. Numerical integration 



The numerical study of eqs. ( |2.5|) -( |2.8D is not as straightforward as one might hope; 
actually we found that, if we start from the UV with initial conditions that reproduce 
the asymptotic behaviour given by eqs. (|3.2| )- (|3.5| ), the numerical integration runs 
almost immediately into divergencies. The reason for this numerical problem will 
become clear at the end of this section, but we have found that starting instead from 
the IR with a solution of the type (|3.11|) and integrating toward the UV poses no 



numerical problem. 

However, this class of IR solutions from which we are starting the integration 
is characterized by more free parameters than the solution (|3.2| )- (|3.5| ) to which we 



would like to connect in the UV. This means that we have to scan the parameter 
space of the IR solutions in order to find those very particular solutions (if any) 
that match to the desired UV behaviour. This can be done as follows. We start by 
studying an IR (p ^> 1) solution of the asymptotic form 

</>{p)=P, £(p)=P, V(P) = ^P + V0, T(p)=0, (4.1) 

keeping r] as the only free parameter. This is a solution of the type fl3.11|) , where 
we have chosen <pi = £i in order to have a confining solution, see eq. (|3.21|) , and we 



have arbitrarily set £i = 1. We also specialize to solutions such that in the IR the 
tachyon goes to zero. 2 Eq. ( |3.13|) then fixes 771 = 1/y/E. The inequalities ( |3.12|) are 



satisfied. The invariance under constant shifts in p allows to fix, e.g., £ — 0. We 
also set 0o — 0; from eq. Q3.19Q we see that a different value of O reflects itself on 



a quantitatively different value of the string tension a, but we expect that this does 



2 We have also studied the case of a tachyon linearly growing in the IR, i.e. t\ 7^ in eq. ( 3.11 ), 
and we have found that in the UV the behaviour of the solution is essentially the same as in the 
case t\ — 0. However, the case t\ = is in a sense more solid because our ignorance on higher 
orders in T of the functions f(T),g(T) becomes irrelevant in the IR. 
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not give rise to qualitative differences in the solution (as we have indeed checked 
numerically, setting 770 = and using £ as a free parameter). 

The results of the numerical integration are as follows. For 770 smaller than a 
critical value r] c , we find that 77 (p) — > —00 at a finite value of p = po, which therefore 
sets a limit to the range of p. The fact that p cannot range from —00 to +00 is to be 
expected. In the usual Dp-brane solutions p is indeed defined only in a semi-infinite 
range which, after a constant shift of p, can always be taken to be < p < 00, and 
p = is the D-brane horizon. However, in our case we find that £(po) is still finite. 
Since £ is related to the energy scale, £ = — 41n(£'/A), this means that in these 
solutions there is no region that we can interpret as an UV limit of the dual theory 
(or, in the gravity language, these solutions have no horizon). Therefore the solutions 
that start in the IR with 770 < rj c are not connected to any of the UV solutions given 



by eqs. (Q-Q. 

If instead 7/0 > 77 c , the situation is reversed and £ — > —00 at a finite value po 
while rj(po) and 0(po) are still finite. This means that at po the energy scale E — ► 00, 
so now the solution reaches the UV regime. However, 0(po) is finite and therefore 
<7ym does not go to zero in the UV. Then, again, this IR solution does not match to 
the desired UV solution. 

The situation is however different for 77 = rj c . In this case <f), £ and 77 all go 
toward —00 at the same value p = p , while T — > — 1; therefore, there are chances of 
matching this solution with the UV solution of eqs. (|3.2| ) -(|3.5| ) . 



It is convenient to introduce the variable y in analogy to eq. (|3.1|), so that the 
UV corresponds to y — > 00. In the analytic computation of sect. @ the UV limit 
was at p = while now it is at p = p , and to have the same convention we must 
make a constant shift in p by — po, so we define y from 

p-po = e~ v . (4.2) 

With repeated runs we have located accurately the value of rj c and the corresponding 
Po, obtaining r] c = 0.4962700(1), po = —0.2076889(1). Locating r\ c with great preci- 
sion is of course necessary if we want to follow the solution deep into the UV regime. 
With a precision O(10 -7 ) on r] c we estimate that we can reliably follow the solution 
in the UV region up to y ~ 7 In 10 = O(10). 

The solution is shown in figs. (|I])-(f|) (solid lines), plotted against y. The func- 
tions £, 77 and T in the UV are very well fitted by the leading terms of eqs. Q3.3|), 



and (|3.5|) , suggesting that we have succeeded in matching the confining IR solution 
to this solution. For the dilaton, however, the situation is more subtle. Fig. @ shows 
e"^/ 2 ~ 1/<7 Y M plotted against — £/4, i.e. against bi(E/A). In the IR the coupling is 
strong and g\ M ~ 1/E 2 , while in the UV regime 1/(? Y m scales linearly with liaE. 

Thus, first of all we see that we indeed succeeded in matching a confining IR 
solution with an asymptotically free UV regime, and now we want to understand 
whether this UV solution is related to eq. (|3.2j). If one compares with the data 
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shown in fig. (|2]) one finds that, numerically, the expansion given by eq. ( |3.2| ) misses 
badly. It is instructive to understand the reason. In the UV region the numerical 
result is very well fitted by 

e -*/a ~A + by + 0(\ny), (4.3) 
and the fit gives b ~ 0.005. Thus, the solution in the UV is well reproduced by 

0~ -2 \n (A + by + (In y)) . (4.4) 
At asymptotically large values of y, this is the same as 

(f) ~ -2 In y - 2 In b + 0(ln y/y) + 0(1 /y) . (4.5) 



This is just the leading behaviour predicted by eq. (|3.2j) , which also predicts b = 
2-15/2 ^ 0.0055, in excellent agreement with the value from the fit. 

However, the expansion of ln(A + by) is valid only if by ^> A, and here A is of 
order one while b ~ 0.0055 is quite small. So, while eq. ( |4.4|) reproduces the data very 
well, its expansion, and therefore eq. ( |3.2| ), is only valid for y ^> 1/6 ~ 181, which is 
well beyond the point where we can push the numerical integration 3 . Then, it should 
be clear that our numerical solution in the UV is indeed nothing but a member of 
the family of solutions given in eqs. (|3.2j)-fl3~5]), and that the analytic expansion for (ft 



given in eq. ( |3.2| ) is only valid for y 3> 181, rather than for y ^> 1. However, to clear 
up any doubt, in the next section we will present an improved analytic solution of 
the equations that is valid for ?/ > 1, rather than only for y 181, that in the region 
y ^> 181 reproduces eq. (|3.2|) , including all the subleading terms written there, and 
that in the region y < 10 where the numerical integration is possible reproduces very 
well the data. An unexpected bonus from this analysis will be that we will also find 
terms ~ exp{— const.?/}, that we will relate in sect. (|6]) to renormalon singularities. 

The physics of the solution is illustrated by figs. (|§)-(0)- Fig. (|5|), as already 
discussed, shows exp{— 0/2}, which is proportional to the inverse coupling l/g^M' 
as a function of \n(E /h). This plot can also be used to define the constant A which 
fixes the energy scale. We define it as the value of E when asymptotic freedom sets 
in, so by definition in fig. (||]) the change of regime takes place at \n(E/A) = 0. On 
the right we have asymptotic freedom, since l/g Y M ~ ^ n (E/A), while on the left side 
of the plot we have strong coupling and confinement. 

Fig. (H) shows the beta functions /3(<7ym) for our solution. It has no zero except 
for the perturbative one at jym = 0, in contrast with the two-loop result, which has 



This also explains why we failed to integrate the solution starting from the UV asymptotics: 



we gave initial conditions that would have reproduced eqs. (|3 . 2| ) — (3.5 ) in a region y ~ 10 where this 
is not a good approximation to the solution. On the other hand, in the numerical integration it is 
very difficult to start from much higher values of y, because one is confronted with exponentially 
small terms in the equations of motion. 
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Figure 1: The result of the numerical integration for T(y) (red solid line) discussed in 
sect. @, together with the analytic results in the UV (black, dashed) and in the IR (blue, 
dotted) discussed in sect. ||[ 



a second zero at l/g Y M = ^V(4tt) 2 , of course well outside the limit of validity of the 
perturbative computation, since there the 't Hooft coupling is Ngy M = (47r) 2 ^> 1. 
The numerical result in the intermediate region matches well the analytic results at 
weak and strong coupling 4 . 

Another interesting quantity is the radius of the 5-sphere (in the Einstein frame, 
since in this frame the UV metric approaches AdS^ x S 5 ) -R(s), which is given by 



eq. (O) 



4) = el"". (4.6) 

If the adjoint scalar fields become massive, the radius of the 5-sphere shrinks to zero. 
We see from fig. (^) that R 2 5 ^ diverges in the IR limit. At the transition between the 
IR and UV region it bounces and finally settles to the constant value 1 / y/2 predicted 
by eqs. ( |3.3|) and ( |3.4| ). We see that it never vanishes, and therefore the 6 adjoint 
scalars remain massless. 



4 With the vertical scale used in the figure, needed to show the matching of the numerical result 
with the IR analytic behaviour, the UV behaviour (3(g) ~ — g 3 cannot be distinguished from the 
horizontal axis, but the numerical result indeed matches with ~ — g 3 in the UV. Note also that, to 



produce this plot, we have set to unit the proportionality constant in eq. ( 2.10 ). Once one has the 
correct proportionality constant, the correct figure is obtained with a rescaling of the units of the 
axes. 
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-0.4^ 



Figure 2: The result of the numerical integration for <fi(y) (red solid line) discussed in 
sect. @, together with the analytic results in the UV (black, dashed) and in the IR (blue, 
dotted) discussed in sect. ||[ 




Figure 3: The result of the numerical integration for £(y) (red solid line) discussed in 
sect. @, together with the analytic results in the UV (black, dashed) and in the IR (blue, 
dotted) discussed in sect. ||| 
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Figure 4: The result of the numerical integration for ij(y) (red solid line) discussed in 
sect. @, together with the analytic results in the UV (black, dashed) and in the IR (blue, 
dotted) discussed in sect. ||| 
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Figure 6: The beta function against gyW- numerical integration (red solid line), and 
IR analytic limit (blue, dotted). The UV analytic limit is (3(g) ~ — g 3 , and cannot be 
distinguished from zero on this vertical scale. 




In (E/A) 



Figure 7: against ln(E/A). 
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5. Analytic solution 



In this section we find an analytic expression for the RG flow which reproduces 
very well the solution from the IR to the UV regions, and even allows to compute 
exponentially small corrections in the UV limit. 

The results of the previous section suggest the following improved UV ansatz 



z -m = A + by + u ( y j 

ay) = -y + ln2 + ((y) 
V 
2 
-1 



v(y) 
T(y) 



+ ln2 +u(y) 
f *(?/)• 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



where u, C, w and t are small, in a sense that will be clarified below. Here 6 = 2~ 15 / 2 ~ 
1/181, and A is a constant of order one, yet to be fixed (A > otherwise the solution 
would not make sense when by < —A). At first, we linearize the equations of motion 
with respect to u, (, uj and t (we will discuss later the non-linear terms), and we get 



2C 



C" + c 

J' + u' 

t" + 1' + It ~ + 



2b 



2u 



A + by 

lb 
A + by 

166 



1 + 



c 



5lu - 2t - 



u 



2t - 



u" + 



1 



2b 



A 



u' + 



by 
b 2 



A + by 
u 

A + by 
u 

— huj — t 

2 A + by 



c 



-u 



c 

6(| - 5u - 2t) . 



(5.5) 
(5.6) 
(5.7) 

(5.8) 



A + by) ' (A + by) 2 ~ A + by 

where the prime is dj dy. We can now take advantage of the fact that b is numerically 
small, and we treat it as a small expansion parameter. Instead A + by in the UV is at 
least of order one. Then, to lowest order in b a particular solution of eqs. ( |5.5| )-( |5~7| ) 
is given simply by 



C(y) = u(y) = -t(y) 



A + by 



1 + 



A + by 



(5.9) 



In fact, in the region where A + by = 0(1), £' is smaller by a factor b compared to 
( while for y » 1/6 we have 1 ^> £ ~ 1/y ^> an d similar for u,t. Substituting 
eq. ( |5.9| ) into eq. ( |5.8| ) we obtain the particular solution 



39 

<y) = -y61n(A + 6y) 



O 



A + by 



(5.10) 



To obtain the most general solution, we first tentatively assume that the terms 
~ C, t and ~ uj (A + by) on the right-hand side of eqs. ( |5.5| )-(|5T8|) can be neglected 
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compared to 1. Then the most general solution is obtained adding to the above 
particular solution the solution of the homogeneous equations 



t" + t' + 2t ~ . 



uj" + J - 2uo ~ , 



u" + u' = . 



(5.11) 



This would give 




+ B T e~ y/2 sin (^-y + «r 
f- B„e~ y , 



(5.12) 
(5.13) 
(5.14) 
(5.15) 



where B^, B n , B?, B u and ay are arbitrary coefficients of the homogeneous solutions. 
Note that the homogeneous equation for the tachyon has two complex solutions 
t = exp{(— 1 ± i\/7)y/2}, and then the solution is an oscillating function, as indeed 
we expected from fig. ([!]). 5 

However, if we now insert back (,uj,t and u/(A + by) on the right-hand side 
of eqs. (|Q| )-(|5"lf) to check whether the approach is self-consistent, we see that the 
tachyon gives a contribution ~ bB T e~ y l 2 sin((-\/7/2)?/ + ar) to the right-hand side of 
so that, for instance, the term e~ 2y in the solution for ( dominates 



HE 



eqs 

over this contribution only if be~ y l 2 e~ 2y or y <C ln(l/6). In the opposite limit 
y ^> ln(l/6) (which with our value of b means i/ > 3) before considering terms 
~ B^e~ 2y we must also include the contribution of the exponentially small term 
coming from t to the right-hand side of eqs. (|5.5p-(|5"T8"D, and we must also go beyond 
the linear approximation: e.g., if we want to get correctly all the terms up to e~ 2y 
we must include the contributions from terms ~ t 2 ,t 3 ,t 4 and u 2 . The correct result 
then turns out to be of the form 



-#v)/2 



39 

A + by - —b \n(A + by) 



O 



A + by 



+ 



bB T e 



-y/2 



sin 



'y/7 



y + a T )[l + 0(b gi (y))} 



(5.16) 



B u + bB 2 g 2 (y)} eT» + O (e" 3 ^ 2 



5 We discarded exponentially growing solutions of eq. ( 5.11 ), as well as the constant solution of 
u" + vf = 0, since the constant part of u(y) has already been taken care of by the constant A, 



eq. (5.1). 
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with gi(y), #2(2/) non-singular and oscillating functions of y. In this work, we are not 
particularly interested in their explicit form (for the comparison with the numerical 
data these terms play a very marginal role), so we computed explicitly only the term 
~ e _y / 2 , at lowest order in b, but it is important conceptually to understand the 
structure of the exponentially small corrections, as we will see in sect. (||). 
A similar structure can be obtained for £,77, T, 



£(v) 



-y + In 2 + 



A + by 



bBn 



>T -y/ 2 ■ [V7 

— —e 91 sin y + a T 

A + by 1 2 y 



V , b 
-£ + ln2 + - — 

2 A + by 



bB 



T{y) 



A + by 



-1 



T e-^ 2 sin 



A + by 



+B T e 




(5.17) 



o( e - y ) 



(5.18) 



0(e- y ) 



(5.19) 



sin I — y + a T 



+ e 



We can now compare with the numerical solution and extract A, Bt, B u , ot from 
a fit. We find A ~ 1.1, B T ~ 0.8, B u ~ — 0.6, a^r — 0.6. This analytic solution is 
shown in figs. (|I|-^) (black dashed line) and we see that it reproduces very well the 
numerical solution already for ?/>l. 

Expanding eqs. ( |5.16| )-( |5.19| ) in the UV we find that they reproduce eqs. 
( |3.5|) , including all subleading terms. The coefficient k is related to A by 



K 



2 A 

39 ln6- — 

b 



585 



In 2 - 256^ A, 



(5.20) 



As we noted above, our solution makes sense only for A > 0, since otherwise \a(A+by) 
becomes imaginary at intermediate values of y. Eq. ( |5.2C| ) then shows that there 
exists a critical value of k, 

k c = -— In 2 (5.21) 



and that our solution (|5.16|) - (|5.19|) can produce, in the deep UV region y 1/6, only 
solutions of the form (|3.2|) - (|3.5|) with k < k c . The solution considered by Klebanov 
and Tseytlin |12| has instead k = and it does not belong to this class. Indeed, 
ref. |12[ finds that in the IR it connects to another conformal fixed point at infinite 
coupling, rather than to a confining solution. 

We see therefore that k plays a crucial role in determining the IR behaviour of 
the RG flow. In the language of Wilson's renormalization group, this means that k 
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is a parameter that determines to which universality class the UV theory belongs, 
and we have found that there are (at least) two universality classes: for k < k c we 
are in the domain of attraction of a confining IR fixed point, while for k > k c we are 
in the domain of attraction of a different fixed point. 

It is also important to understand how the result changes if we change the form 
of the function f(T). We find that, as long as /(T) has a minimum at a finite 
value Tq, the qualitative behaviour of the solution does not change. However, it is 
especially important to known whether the exponential terms ~ e~ y ^ 2 are affected 
by the form of /(T). Linearizing eq. ( |2.8| ) for f(T) generic, and writing T = To + 1, 
the left-hand side of eq. (0) becomes t" +f + ct, with c = /"(T )/(2/ 2 (T )); for 
f{T) = 1 + T + T 2 /2, c = 2. The homogeneous equation now has the solutions 
t = e ty with 



7 = --(l±zV4c^T). (5.22) 

We see that, as long as c > 1/4, Re 7 is unchanged and we still have solutions of 
the form e~ y l 2 sin(z/y), but the frequency v depends on the explicit form of f(T). 
Therefore the value that we found, v = a/7/2, does not have a special significance, 
being related to the unknown form of f(T), while the factor e~ y l 2 is universal as long 
as f(T) obeys 

2/"(T ) > / 2 (T ) . (5.23) 

The function g(T) instead does not enter the homogeneous equation, so from this 
point of view its form is irrelevant. 

Finally, it is straightforward to compute analytically the corrections to the solu- 
tion in the IR regime. The corrections are exponentially small in p, 

0=0 1 p + 0o-E^ )ex P{-7 l W p} ) 

i=l 

£= £iP + £o -£ Cf exp{- 7 p ] p} , (5.24) 

i=l 

V = V1P + V0 - £ Q exp{-7^ ) p} , 

i=l 
i=l 

We have computed the coefficients Ci, 7, for i = 1,2, 3. We get, for 0i = £1 = Vo^i, 

# = (V5 + 3)0! 7? = 2(75 + 1)0, 7 f = ^f^i 
7 P = 20! = 40! T f = (V5 + 3)0! 

7^ = # = ^01 7^ = (V5 + 3)0! 

7P = 20! 7 f } = (y/5 + l)0i 7 f ) = 40i 
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(jW _ \ e |0o-|?o-5r;o (jW \ </>o-3£o-10»7o 

1 " 16(75 + 3)20? 2 16(V5 + 1) 4 0! 

rtf) _ 125 c i<t>o-Uo-9r,a _ 1 c -2g 

3 32(9 + 3V5) 2 0? 1 20? 

/7«) = _JL p -4?o r (0 _ 1 f 2 M p^-^o-Srjo 

2 " 160f ° 3 ^ (v/5 + 3)20f ^(^5 + 1)2 • 1 
Gl "20f e ° 2 "40f 6 
63 " 16(^5 + « 6 Cl " 202 e 

(T) = 1 e ^-|fa-5np £(21 = 3 -4&, 

2 (75 + l) 2 0f 3 1601 




This analytic solution is shown in figs. as a blue dotted line, and we see that 

it reproduces very well the numerical solution in the IR region. Having included the 
the first three subleading exponential, we see that the IR analytic solution has an 
overlap with the analytic UV solution, so that the two expansions together describe 
analytically the solution in the whole region. 



6. The renormalon singularities 

A striking feature of the UV solution is the appearance of the exponentially small 
terms ~ e~ y ^ 2 , e~ y , etc.; in the UV, at leading order, y ~ — £ = AlnE, so (extracting 
an overall proportionality factor 46) eq. ( |5.16 ) gives 

1 Inf -?lnlnf + ^F l{ E)^ + B 2 F 2 (E)^ + C^ + 0(±-J (6.1) 



g\ M {E) '"A 8 A ' 4 ~ iy ~ J E 2 ' "'"^-'£4 1 ~ #4 ■ ~^ E e 
with 

Fl (E)=sm((3 T \nj + a T \ [1 + 0(6)], (6.2) 

and C = B u j (46) (however, we expect that the numerical values of these coefficients 
will be affected by our uncertainties, as it happens to the factor 39/8). F 2 (E) can 
be determined adding the terms ~ t 2 to eqs. ( |5.5| ) and ( |5.8| ). The positive constant 
Pt is given by 



2/"(T ) A 1 ' 1 



(6.3) 



/ 2 (T ) 

and depends on the explicit form of f(T), so its numerical value is not under control. 
The exponential correction to the relation y ~ — £ gives an 0(6) contribution to 

G3)- 

First of all, we read from eq. ( |6.1| ) that there are power corrections, and they are 
exactly of the form that is expected from the effect of renormalons. 
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Let us recall that, given a divergent series 



G(g 2 YM ) = E G n g 2 Y n u , (6-4) 

n=l 

with \G n \ bounded by r n n\ at large n, one defines its Borel transform B(z) as 

oo /~t 

B{z) = Y J ^z\ (6.5) 

n=0 U - 

inside its convergence radius, \zr\ < 1, and by analytic continuation elsewhere. Then 

G(g 2 YM )= / dze-'^B(z) (6.6) 
Jo 

(where the integration contour is on the positive real axis) can provide a resummation 
of the perturbative series, depending on the structure of the singularities of B(z) in 
the complex plane (Borel plane), and on the convergence properties of the integral at 
infinity. For asymptotically free non abelian gauge theories a number of singularities 



in the Borel plane have been identified as follows [^, |30|, [31], |32|]. IR renormalons 
give singularities on the positive semiaxis, at 

z=%-k, k = 2,3,... (6.7) 

Po 

(po > for asymptotic freedom; in our case fi = |iV/(87r 2 ), see eq. (|3.9|)). Since 
these singularities are on the positive semiaxis, we must supplement eq. (|6.6|) with a 
prescription for dealing with them, and this produces non-perturbative contributions 
of order 

2k 1 

exp{ -^ } ~^' (6 ' 8) 

times regular functions of the energy. Since the first IR renormalon is at k = 2, it 
gives a contribution ~ 1/E 4 . UV renormalons are instead on the negative semiaxis, 
at 

z = --k, fc = l,2,3,... (6.9) 

Po 

Even if these singularities are not on the integration contour, they limit the conver- 
gence radius of the expansion for B(z) and are responsible for effects of the form 

exp{--^} ~ — I (6.10) 
3ym & 

However, now k starts from 1, so that the leading effect is 1/E 2 . There are fur- 
ther singularities on the positive semiaxis due to instantons, but their effect is 
of order exp(— 87r 2 /g Y M)5 which for SU(N) with 6 adjoint scalars is ~ (l/E) 8N ^ 3 
(or ~ (1/E) 16N / 3 if the first non- vanishing contribution comes from an instanton- 
antiinstanton pair) and therefore is quite suppressed for large N. 
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Comparing our result ( |6.1| ) with the above situation, we see that the 1/E 2 term 
that we have found matches with the contribution of the first UV renormalon. The 
fact that we get the correct location of the renormalon singularity in the Borel plane 
is non-trivial, and can be traced back to the fact that the linearization of eq. ( |2.S| ) 
gave eq. ( |5.7| ), which has the homogeneous solution t(y) = e lv with Re 7 = —1/2, 
as long as eq. ( |5.23| ) holds. From the point of view of the gravity computation, we 



could have obtained 1/E a with a a priori any real number, not necessarily integer 
nor rational. We therefore regard this agreement as a successful and non trivial test 
of the conjectured duality between type OB theory and non supersymmetric YM 
theory. 

The term B^F2(E)/E A comes from the iteration of the term BtFi(E)/E 2 , as 
we see from the coefficient B^, so it clearly corresponds to the UV renormalon with 
k = 2. Iterating further the contribution of the tachyon in the solution of the 
equations of motion, we find all the UV renormalons at their correct locations in the 
Borel plane. 

The term C/E A instead has a prefactor which is just equal to the constant C and 
is not related to Fi(E) or B?, so it has a different physical origin. Thus, we identify 
it with the k = 2 IR renormalon. Again, it is very remarkable that its position in 
the Borel plane matches exactly the position of the first (k = 2) IR renormalon; 
furthermore, we correctly find no term corresponding to an IR renormalon with 
k = 1, which corresponds to the absence of gauge-invariant operators of dimension 
2 in the OPE of current- current correlation functions. 

Finally, even the energy dependence of the prefactors is consistent with what 
is known about renormalons from the analysis of Feynman graphs. In fact, in the 
approximation in which only a single chain of bubbles is considered, it is known that 
all UV renormalons are double poles, while the k = 2 IR renormalon is a single pole, 



and all others IR renormalons are again double poles ||32| . For a single pole 



B(z) ~ — r, (6.11) 

Z A) 

and then eq. ( |6.6| ) gives a contribution ~ exp{— 4/(/?o5'ym)} ~ 1/E 4 " without any 
further energy dependence, in full agreement with our interpretation of C/E A as the 
contribution of the k = 2 IR renormalon. For a double pole on the positive semiaxis, 
instead, the non-perturbative contribution is of order 

/•oo 9 I 

G(g 2 YM ) ~ / dze'^ru 



[z — z) 2 



1 r°° 12 1 1-/2 
-> -5- / dze~ z ^M ~ e -*/9 YM ( 6-12 ) 

9ym j o - z) 9ym 

where in the second line we have integrated by parts and kept only the non-analytic 
term. This provides a InE term in the prefactor. Similarly, there are logarithmic 
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Figure 8: The function F%{E). 



terms from the double poles on the negative semiaxis |32|, |33| . Thus in the single 
chain approximation the first UV renormalon is of order 

~-^> (6-13) 

and in our approach is reproduced by the first term of the expansion of eq. ( |6.2j ) 
in powers of (3t (which suggests that (3t is proportional to /3o). Higher order terms 
are expected to come from graphs with multiple chains insertions, and are not well 
understood in QCD fl32| . In this sense, eq. (|6.2|) is also a remarkable prediction of 
the dual string theory, that would be very interesting to check from the YM side. 

At low energies F\(E) oscillates faster and faster, see (fig. |8|). This is physically 
very reasonable, since these oscillations signal the transition from asymptotic freedom 
to the confining regime. Thus, the form of Fi(E) seems to capture both the correct 
single chain result and a physically correct behaviour when the IR is approached. 

We conclude this section with a general comment. The rule concerning the posi- 
tion of renormalons in field theory is that we find on the positive semiaxis the renor- 
malons corresponding to strongly coupled physics. For asymptotically free gauge 
theories IR renormalons are at z > 0, while UV renormalons are at z < 0. In QED 
or A0 4 the coupling becomes strong in the UV and correspondingly the UV renor- 
malons are at z > and IR renormalons are at z < 0. When a singularity is at z > 
we must give a prescription for going around it; if we understand the physics asso- 
ciated to this singularity, we can give a physically motivated prescription; then we 
are left with a non-perturbative contribution, but no ambiguity on the resummation 
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of the perturbative series. This is the case of instantons . If however we do not 
understand the physics behind the singularity, we do not know what prescription to 
use. Different prescriptions (e.g. going around a pole from above or from below) 
give answers that differ by terms of order exp(— 2k/(3og\ u ), so renormalons on the 
positive semiaxis are usually taken as a source of an uncertainty of this order in the 
resummation of the perturbative series 6 . 

In our case, however, the coefficients Bt, C in eq. (|6.1| ) are fixed from the com- 
parison with the numerical solution. Once we specify the IR solution, the values of 
Bt, C, etc. follow, and the renormalon contribution, rather than an uncertainty on 
the perturbative series, becomes a well defined non-perturbative contribution. The 
fact that IR physics fixes the renormalon contributions confirms the old conjecture 
of 't Hooft that renormalon singularities are related to the quark confinement 
mechanism. It is interesting to note that instead Bt, C cannot be fixed, even in 
principle, from the UV side, since there they just appear as integration constants of 
the homogeneous equations (|5.11|) . 



7. a' and string loop corrections 

We next consider the effect of a' corrections on the solutions. For the type theory, 
these have been discussed in ref. [|IT|, where it is found that they have the same 
structure as in type II theories; in particular, world-sheet supersymmetry implies 
the vanishing of ~ a'(Riemann) 2 corrections, and the first non-vanishing correction 
is ~ a /3 (Riemann) 4 . It is well known that a' corrections suffer from a certain degree 
of ambiguity, because not all the coefficients of the possible operators are fixed by the 
comparison with the string amplitudes; a related ambiguity is the fact that one can 
perform field redefinitions that mix different orders in a', e.g. g^ u — > g^ u + a'R^ u . If 
one has an exactly conformal background these field redefinitions do not matter, but 
of course if one works at a finite order in the a' expansion they make a difference. 

An important observation |12] is that we can choose the a! correction such that 
they depend only on the Weyl tensor, rather than on the Riemann tensor, and then 
they have a much milder behaviour when we approach AdS§. To be quantitative, 
consider eq. Q2.5Q , which, including the first non- vanishing a' correction, has the 
general form [|12| 

4> + l -g{T)e^-^ = ce^" 5 " (Weyl) 4 , (7.1) 

where c is a constant, (Weyl) 4 denotes the appropriate contractions of four Weyl 
tensors Cmnrs, and again we have set a' = 1. Writing (p = e 2y ((p" + </>'), where as 



6 Actually, in QCD even UV renormalons can give rise to uncertainties on the perturbative series 
which from the Feynman graphs point of view are uncalculable, despite the fact that the UV physics 



is well understood |34 
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before the dot is d/dp and the prime is d/dy, with p = e y , we have 

0" + 0' + I^(T)e^- 5 "-^ = ce^" 5 ^ (Weyl) 4 . (7.2) 



We see that an ansatz of the form ( p.2|) -( p3|) , or its improved form ( |5.1|) -(5.4), is still 



a solution, but the numerical coefficients in front of the factors In 2 change. In fact, 
consider the ansatz (p = — 2 hay + a\ In 2 + O (In y/y), £ = —y + a 2 In 2 + 0(l/y), r\ = 
—y/2 + a 3 ln2 + 0(l/y),T = —1 + 0(l/y). We have checked that on this ansatz all 
components of the Weyl tensor Cmnrs are 0(l/y), so that W2\ 



(Weyl) 4 ~i_. (7.3) 
Substituting into eq. ( |7T2] ) we get (setting again g{T) = T 2 /2) 

_2 + 1 e (^-5 a3 )ln2 0( 1 } = OOTSt e(z 3^_ 5a3)ln2 1 

and similarly for the equations for £, r/, T. Thus, the leading y dependence cancels, 
and we can adjust the coefficients 01,02,03 so that the equations are satisfied. We 
see that the numerical values of the coefficients ai, a 2 , a 3 are affected by the inclusion 
of the (Weyl) 4 term. The logarithmic behaviour in the UV is not spoiled, but the 



numerical values of the beta function coefficients change [|T^ . 

We now want to see whether the terms related to renormalons are affected. If 
we use again the improved ansatz ( [5.1| )-( |5~4D , where now b and the coefficients of 
In 2 in £,77 are modified according to the discussion above, we see that, in the limit 
by ^ 1, where b/(A + by) ~ 1/y, the net effect of the addition of the (Weyl) 4 term is 
to change the numerical factors inside the square brackets on the right-hand side of 
eqs. (|5.5|)-(|5T8l). However, the associated homogeneous equations are still given by 
eq. ( |5.11|) , and are unaffected by the inclusion of the (Weyl) 4 term. Therefore, the 
position of the renormalon singularities in the Borel plane is unaffected by these a' 3 
corrections. It is clear that the argument goes through at all orders in a' as long as 
the corrections can be written in terms of the Weyl tensor. As for the string loop 
corrections, in the UV they are small by definition, since e*<l. Thus, it appears 
that the result on the position of renormalon singularities is quite robust. 

In the IR the situation is more complicated, first of all because unavoidably 
when oym becomes strong also the string loop corrections become important. This is 
different from the standard AdS / CFT correspondence, where the dilaton is constant 
and g s = can be taken to be small everywhere. 



The role of a' corrections in the IR has been discussed in ref. [14|, and it has been 



found that they depend crucially on the numerical relations between the coefficients 



£1, 01, 771 that parametrize the IR solution (|3.11|). With the choice £1 = <pi that we 



used in sect. f|, the corrections are indeed large in the IR, while, for instance, for 
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Figure 9: e ^l 2 and the tachyon in the case (pi = 3, £1 = 1. 
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Figure 10: R 2 5 ^ in the case 4>\ = 3, £i = 1. 

0i = 3£i,£i = we still have a confining solution but now the a' corrections in the 
IR are exponentially small in p. Therefore we have repeated the numerical analysis 
discussed in sect. |], for the case 4>i = 3£i, £i = 1, t\ = and, from eq. ( |3.13| ), r\\ = 1. 
We have found that the lowest order solution for the dilaton and the tachyon shows 
no qualitative difference, see fig. (fj|). 

The functions r], £ are also qualitatively very similar to the previous case; the 
main difference is that, since now ^i/2 — 771 < 0, the radius of the 5-sphere R 2 5 ^ = ei~ n 
goes to zero in the IR, see fig. fllCf) (using eq. ( |3.13| ), we see that this happens in 
general if (f)\ > (3/2)£f). This suggests that the six adjoint scalars become massive, 
but to clarify this point we need a better understanding of string loop corrections in 
the IR limit. 
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